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ABSTRACT 

Broad Fe Ka emission lines have recently been reported in a number of Galactic black holes. Such lines 
are useful accretion flow diagnostics because they may be produced at the inner accretion disk and shaped by 



relativistic effects, but in general they have only been observed at luminosities of Lx 



10 



37-38 



erg s in soft 



X-rays. The Galactic microquasar V4641 Sgr — widely known for its 12.2 Crab (1.5-12 keV) outburst in 1999 
September — displayed low-level activity in 1999 March. BeppoSAX observed the source in this state and Fe Ka 
line emission was found (in 't Zand et al. 2000). In re-analyzing these data, we find strong evidence that the 
Fe Ka line profile is broadened. For the most likely values of the source distance and black hole mass measured 
by Orosz et al. (2001), our fits to the total spectrum indicate that the source was observed at a luminosity of 



L x = l.9!k% x 10 36 erg s" 1 (2-10 keV), or L x /L 



Eddington 



= l-8±!i o x 10 3 . Advection-dominated accretion flow 



(ADAF) models predict a radially-recessed disk in this regime. In contrast, fits to the observed Fe Ka emission 
line profile with a relativistic line model constrain the inner disk to be consistent with the marginally stable circular 
orbit of a Schwarzschild black hole. 



(60° < i < 71°), and a large distance (7.40 < d < 12.31 kpc, 
Orosz et al. 2001). Hjellming et al. (2000) monitored the 
outburst in radio and found jets expanding at apparently su- 
perluminal velocities; V4641 Sgr is therefore termed a "mi- 
croquasar." Significant low-luminosity activity was observed 
by BeppoSAX and by RXTE Galactic bulge scans (Swank & 
Markwardt 2001) earlier in 1999, however, well before the 
better-known outburst. In 't Zand et al. (2000) reported on a 
pointed BeppoSAX observation made on 1999 March 13. Aline 
complex was noted, and attributed to Fe XXV and Fe XXVI 
lines with widths narrower than the BeppoSAX resolution. The 
best-fit time-averaged 2-10 keV flux measured on March 13 
was found to be 2.7 x 10~ 10 erg cm" 2 s , corresponding to 
1.8 x 10 36 erg s" 1 < L x < 4.9 x 10 36 erg s _1 . Motivated by the 
constraints on the accretion flow geometry made through the 
study of broad lines in other sources at higher luminosities, we 
re-visit this intermediate-luminosity spectrum of V4641 Sgr. 

2. OBSERVATION AND DATA REDUCTION 

BeppoSAX observed V4641 Sgr on 1999 March 13.22. 
The soft-X-ray-imaging narrow-field instruments (NFI) aboard 
BeppoSAX are called the Low-Energy and Medium-Energy 
Concentrator Spectrometer (LECS and MECS, respectively). 
The net exposure time for the LECS was 10.2 ks; the net MECS 
exposure was 27.7 ks. The MECS has an energy resolution of 
8% at 6 keV. 

The Phoswich Detector System (PDS; 12-300 keV) ob- 
served V4641 Sgr simultaneously with the LECS and MECS. 
The statistical quality of the high energy spectrum obtained is 
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1. INTRODUCTION 

Recently, studies of the Fe Ka line profiles in Galactic black 
holes and black hole candidates have revealed strong, skewed 
lines, similar to those that are observed in some Seyfert galax- 
ies (see, e.g., Fabian et al. 2000; Weaver, Gelbord, & Yaqoob 
2001). Such lines are likely produced by the irradiation of an 
accretion disk which extends close to the marginally stable orbit 
around the black hole. The exact shape of such lines can indi- 
cate the inner extent of the disk (and thereby constrain black 
hole spin), its ionization state, and the nature of the interaction 
between the disk and the corona that irradiates its surface. 

The broad Fe Ka lines recently reported have all been ob- 
served in bright persistent systems (Cygnus X-l: Miller et 
al. 2002a, GRS 1915+105: Martocchia et al. 2002) or in 
bright states of transient systems (XTE J 1650-500: Miller et 
al. 2002b, SAX J171 1.6-3808: in 't Zand et al. 2002). These 
lines provide important new constraints on the black hole spin 
and inner accretion flows of these systems. While most mod- 
els for the accretion flow geometry in Galactic black holes pre- 
dict that the disk extends near to the marginally stable orbit in 
bright states (Lx — 10 37-38 erg/s, or higher), the nature of lower- 
luminosity states may be different. In particular, ADAF models 
for transient systems (e.g., Esin, McClintock, & Narayan 1997) 
predict that the disk may be recessed to several tens or hundreds 
of gravitational radii (R g = GMbh/c 1 ). 

V4641 Sgr is usually identified with its remarkable 12.2 Crab 
outburst lasting only a few hours on 1999 September 15 (Smith 
et al. 1999). Dynamical constraints strongly indicate a black 
hole primary (5.49 < M BH < 8.14 M ), a high inclination 



V4641 Sgr 



relatively poor, and does not allow for tighter constraints on the 
continuum. In particular, reflection models — which describe 
a geometry in which the Fe Ka emission line and a "Compton 
hump" result from the irradiation of the disk — are not better- 
constrained by considering the PDS spectrum. Our analysis is 
therefore limited to the time-averaged NFI spectra. 

All aspects of the source and background spectral extraction, 
instrument response function generation, and filtering are ex- 
actly the same as reported in in 't Zand et al. (2000). We 
considered the LECS spectrum between 0.4 and 4.0 keV, the 
MECS spectrum between 2.0 and 10.0 keV (again, similar to 
in 't Zand et al. 2000). The spectra were fit using XSPEC ver- 
sion 11.1 (Arnaud 1996). All errors reported in this work are 
90% confidence errors. The LECS and MECS spectra were fit 
jointly. An overall normalizing constant was allowed to float 
between the spectra to account for differences in the instrumen- 
tal flux calibration. 



3. ANALYSIS AND RESULTS 

Preliminary fits suggested that the absorbing column and the 
nature of the soft component could not be constrained inde- 
pendently. We therefore fixed the equivalent neutral hydro- 
gen column density at the weighted average measured along 
this line of sight by Dickey & Lockman (1990): N H =2.3 x 
10 21 atoms cm -2 (assuming the "phabs" model in XSPEC). 

We attempted to describe the continuum spectra with a num- 
ber of models, including the "bulk motion Comptonization" 
model (Shrader & Titarchuk 1999), a model consisting of a 
multicolor disk black-body (MCD; Mitsuda et al. 1984) and 
power-law components, and a model consisting of MCD and 
"comptt" (Titarchuk 1994) components (as per in 't Zand et al. 
2000). None of these models gives an acceptable fit to the data, 
yielding y?/v = 6.22, 5.83, and 5.96, respectively (v is the 
number of degrees of freedom; v = 16 for the first two models, 
v = 74 for the last). These continuum models all fail to account 
for a strong emission feature in the Fe Ka line region. 

In Figure 1 , the line profile is shown as a ratio of the spectrum 
to these continuum models. Following Iwasawa et al. (1999), 
the ratio is formed by ignoring the 4-7 keV range when fitting 
the spectra. It is apparent that the line may not be intrinsically 
narrow — the profile revealed in the data/model ratio is similar 
to those expected for a line produced by irradiation of the inner 
accretion disk around a black hole (Fabian et al. 1989, Laor 
1991). In a previous examination of this spectrum, in 't Zand 
et al. (2000) modeled the Fe Ka line region with two Gaus- 
sians (at 6.68 keV and 6.97 keV, as per Fe XXV and Fe XXVI) 
— each with zero widths — assuming that the intrinsic width 
of each line is below the instrumental resolution. The bottom 
panel in Figure 1 shows the data/model ratio obtained with this 
narrow line model; clearly, flux between 5.0-6.4 keV — per- 
haps associated with an intrinsically broad line — is not ac- 
counted for by this model. 

Adopting the MCD plus power-law model for the continuum, 
we next investigated the nature of the Fe Ka line region in de- 
tail. For a line model consisting of two zero-width Gaussians, 
X 2 /V = 1 .721 , (v = 74). However, for fits with the Laor 
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FIG. 1. — Data/model ratios for the MECS spectra and a number 
of viable of continuum models. In black (and with diamonds), the ra- 
tio for an MCD plus power-law continuum. In red, the ratio for an 
MCD plus "comptt" continuum. In blue, the ratio for the "bulk mo- 
tion Comptonization" model. Above: The 4-7 keV range was ignored 
in fitting the spectrum. The red wing of the line is clearly evident 
in the ratio of the data to each continuum model. Below, two zero- 
width Gaussian lines added at 6.68 keV and 6.97 keV were added to 
the models (corresponding to narrow Fe XXV and Fe XXVI emission 
lines). The red wing between ~ 5-6.4 keV is not accounted for by 
these Gaussians. 

line model, y 2 l v = L230,(z/ = 71). The improvement in the 
fit statistic is significant at the 4.5er level of confidence. This 
continuum plus line model is our best-fit spectral model, and 
it is shown in Figure 2. The inclusion of a smeared Fe K edge 
("smedge," Ebisawa et al. 1994) does not improve the fit signif- 
icantly. When a smeared edge is included, the edge parameters 
are not well-constrained: E e d ge = 9.1^' 2 keV and r = l.b^ 1 ^ 
Allowing a standard narrow edge component to float in the 
Fe Ka region does not give an improvement in the fit statistic. 

The Laor line energy was only permitted to vary within the 
range 6.40 keV < E^,,,. < 6.97 keV (Fe I - Fe XXVI). The 
best-fit energy is measured to be Eia 0r = 6.40 +0,3 keV. The line 
emitting region is relatively well constrained: Rj„ = 8^ R g and 
Rou, = 30^ R g (R g = GMbh/c 1 ). As is evident fr om the line 
profile shown in Figure 1, the line equivalent width is very high: 
W Ka = 870^° eV (9.8!^ x 10" 12 erg s" 1 , 0.4-10 keV). The in- 
clination is measured to be / = 43° ± 15°. This inclination is 
lower than the optically-determined range; however, Orosz et 
al. (2001) note that the inclination implied by the jet is very 
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FIG. 2. — Above: The unfolded 0.4-10.0 keV NFI spectrum and best- 
fit model (the spectra shown are normalized to have the same flux). In 
red, the total spectral model is shown. The Laor line model is shown 
in blue. Below: The data/model ratio for best-fit model shown above. 

low (i < 6°). In this context, the intermediate value we mea- 
sure is acceptable. The inner disk emissivity profile (J(r) ~ r~@) 
measured with this model is consistent with that expected for a 
standard accretion disk (J3 = 3). 

With the Laor model for the Fe Ka line region, we find that 
the MCD inner disk color temperature is: kT = 1 . 14 ± 0.03 ke V 
with a normalization of Nmcd = 4.5 ± 0.5; this component con- 
tributes an unabsorbed flux of 1.5 ±0.2 x 10~ 10 erg cm" 2 s" 1 
in the 0.4-10 keV band. The power-law index is measured 
to be T = 1.3 ±0.1 with a normalization of N p i = 1.1 ±0.3 x 
10~ 2 (photons keV -1 cm" 2 s" 1 atl keV); this component con- 
tributes an unabsorbed flux of 1.1 ±0.2 x 10~ 10 erg cm" 2 s -1 in 
the 0.4-10 keV band. 

For the most likely values of the source distance and black 
hole mass measured by Orosz et al. (2001), we measure 
a total source luminosity of L2-10 = 1-9^q'q x 10 36 erg s" 1 , 
or Lo.4-10 = 3.2+};] x 10 36 erg s" 1 (0.4-10.0 keV). As a frac- 
tion of the Eddington luminosity, these values correspond to 
L2-10/W = 1.838 x 10 " 3 and LoA-w/L Edd . = 3.0t\i x 10" 3 , 
or riiEdd. ^ 3.0+[;f x 10" 2 (defining m Edd . = 10 L ho i ometric /c 2 ; 
equivalent to assuming an efficiency of 10%). Extending these 
values to a broader range depends slightly on the spectral model 
used, and the data at high energies is not very constraining. For 



the model adopted by in 't Zand et al. (2000) and our best-fit 
model, the 0.4-120 keV luminosity is approximately double the 
0.4-10 keV luminosity. 

Although fits with three continuum models commonly ap- 
plied to Galactic black hole spectra indicate a broad Fe Ka 
emission line, we have found two sets of fits that indicate 
instead a narrow line. The first set of fits invoke a partial 
covering model. For example, modeling absorption in the 
ISM with "phabs" (N h .ism = 4.8!^ x 10 21 cm" 2 ) and allow- 
ing two partial covering models ("pcfabs" in XSPEC; Nh,\ = 
4.5 ± 1 .5 x 10 22 cm" 2 and ft = 0.35 ± 0.10, N H ,2 = 2.7 ± 0.3 x 
10 21 cm" 2 and fi = 0.66 ± 0.05) and a simple power-law contin- 
uum (T = 2.7 ±0.1) gives an acceptable fit ( X 2 /v = 1 . 10, v = 70) 
with only a narrow Gaussian line (E = 6.89 ±0.04 keV, FWHM 
= O.l+o } keV). Partial covering models are commonly fit to 
neutron star "accretion disk corona" sources which display 
strong dipping behavior in their X-ray lightcurves (see, e.g., 
Morley et al. 1999). Dips have not been reported in V4641 Sgr. 
It is unlikely that partial covering models are justified for this 
source. 

The second set of fits involves the "compps" Comptonization 
model of Poutanen & Svensson (1996). For the models dis- 
cussed above that utilized separate soft and hard phenomeno- 
logical fit components, there is some worry that the red tail of 
the broad line is an artifact of the region where the soft and 
hard fit components cross. We have used compps to model a 
disk blackbody spectrum Comptonized by a slab corona (with 
unity covering fraction), and reflected off of a neutral, cold 
medium with Fe abundance relative to hydrogen that is three 
times solar. An w 40keV corona with r es « 0.9 fits the contin- 
uum well, whereas the line region is well-described by a narrow 
6.9 keV line with w 250 eV equivalent width. However, an ex- 
tremely large reflection fraction of R w 3 is then required. A 
very large neutral Fe edge from reflection coupled with a nar- 
row Fe XXVI line or K/3 line from more neutral Fe would be 
fairly unusual. We note, however, that a very similar compps 
model with a somewhat lower reflection fraction (R w 1 .5), does 
allow for a relativistic line with rest energy 6.9 keV, emissivity 
index fss -2.2, equivalent width 360 eV, and an inner emission 
radius of 6GM/c 2 . The lack of reliable data above lOkeV pre- 
vents us from exploring further these models requiring large 
reflection. 

It is difficult to identify the "state" in which this observa- 
tion of V4641 Sgr occurred (for a recent discussion of Galac- 
tic black hole states, see Homan et al. 2001; for a review see 
Done 2002). In 't Zand et al. (2000) report no measurable 
narrow or broad-band features in the power density spectrum 
(PDS) of the MECS light curve, consistent with the "high/soft" 
state. This state is defined by a soft, thermal, disk-dominated 
spectrum, with a weak, soft (T ~ 3) power-law; we measure 
T = 1.3 ±0.1. The intermediate luminosity and relatively hot 
accretion disk (kT = 1.14 ±0.03 keV, representing 58% of the 
0.4-10 keV flux) are inconsistent with a "low/hard" state iden- 
tification. The lack of broad-band noise in the PDS and the 
extremely hard power-law index are inconsistent with both the 
"very high" and "intermediate" states, though the limits on the 
noise are not very constraining. We suggest that V4641 Sgr 
was observed in a kind of "intermediate" state, which may dif- 
fer from similar states in other Galactic black hole systems. 



4. DISCUSSION 
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We have re-analyzed the BeppoSAX NFI spectra of the 
Galactic microquasar V4641 Sgr obtained on 1999 March 13. 
Applying the Laor model for line emission from the inner part 
of an accretion disk around a black hole, we find that a broad, 
relativistically-shaped Fe Ka line is preferred over a sum of 
narrow helium-like and hydrogenic Fe Ka lines at the 4.5a 
level of confidence. A line model — whether a complex of 
narrow lines or a broadened line — is strongly required by the 
data. Within the context of the Laor model, the narrow emis- 
sion observed near 6.8 keV (see Figure 1) is the blue wing of 
a relativistically-shaped line from more neutral Fe species in- 
cluding Fe I Ka at 6.40 keV (see Figure 2, Section 3). The line 
is among the strongest yet observed in a Galactic black hole 
system: Wx a = 870^'q° eV. A spinning black hole (implied if 
Ri„ < 6 R g ) is allowed but not required by the line model; our 
best-fit value is R in = 8^ R g . Although the line may be broad- 
ened partially by Compton processes, it is likely that relativistic 
shifts dominate as the profile is strongly asymmetric. 

In most cases, observations of Galactic black holes at Lx — 
10 36 erg s" 1 do not achieve the sensitivity required to clearly 
detect Fe Ka emission lines or fast quasi-periodic oscillations 
(few x 100 Hz; plausibly associated with the Keplerian fre- 
quency near the marginally stable orbit). The Fe Ka emission 
line we have observed in V4641 Sgr, then, represents a rare 
opportunity to test accretion flow models. 

Esin, McClintock, & Narayan (1997) predict that at riiEdd. < 
0.08, the inner accretion flow in Galactic black holes will be 
radiatively inefficient and the disk will be recessed to /?;„ > 
200 R g , The line we have observed, however, is consistent with 
being produced in a disk with /?,„ = 8*3 R g . In the 0.4-10.0 keV 
band, we measure riiEdd. = 3.0^ 3 x 10~ 2 (assuming the most 
likely values for the black hole mass and distance measured by 
Orosz et al. 2001). This value approximately doubles on the 
0.4-120 keV range; the value of riiEdd. given by a bolometric 
luminosity might be a factor of a few higher. We also note that 
taking extremes for the black hole mass, distance, and observed 
flux can yield values of riiEdd. that are significantly lower than 
than the range we quote, and also above the critical value of 
riiEdd. = 0.08. It is clear, however, that in the regime where an 
ADAF is expected to begin to take hold, the observed Fe Ka 
line profile strongly suggests a disk which extends close to the 



marginally stable orbit. 

The Fe Ka line we have measured is not the strongest yet 
observed from V4641 Sgr. RXTE/PCA spectra obtained during 
the decline of the 12.2 Crab flare (L x ~ 3-4 x 10 39 erg s" 1 ) 
in 1999 September, reveal an Fe Ka line at ~ 6.70 keV with 
Wko ~ 2.4 keV (Revnivtsev et al. 2002). The resolution of the 
PCA is not well-suited to reveal a line profile of the kind we 
find in the 1999 March NFI spectra. In part, the large equiva- 
lent width we have measured may be explained by an ionized 
disk (Ballantyne, Ross, & Fabian 2002). We also speculate that 
Fe may be over-abundant (relative to solar) by a factor of a few 
in this source. Orosz et al. (2001) report that N and Ti may be 
over-abundant by a factor of 10, Mg by a factor of 7, O by a 
factor of 3, and Ca by a factor of 2. 

The Fe Ka emission line observed in an XMM-Newton spec- 
trum of the Galactic black hole candidate XTE J 1650-500 ex- 
tends to ^4 ke V and implies a spinning black hole (Miller et al. 
2002b). The line profile seen in the NFI spectra of V4641 Sgr 
is not as strongly skewed (perhaps extending to ^5 keV; see 
Figure 1) as the profile observed in XTE J 1650-500, likely 
indicating that the line is shaped less by strong gravitational 
effects. In an NFI observation of the Galactic microquasar 
GRS 1915+105, Martocchia et al. (2002) find a relativistically- 
shaped Fe Ka line consistent with neutral ion species, which 
does not require a spinning black hole (based on the implied 
inner extent of the accretion disk). Qualitatively, the line pro- 
file observed in GRS 1915+105 may be the most similar to that 
in V4641 Sgr. RXTE/PCA spectra of the Galactic black hole 
candidate GX 339-4 at Lx ~ few x 10 36 erg s" 1 also reveal ev- 
idence for broad, relativistically-shaped line profiles (Nowak, 
Wilms, & Dove 2002) that render an ADAF geometry unlikely. 
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